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Aba^aract 

A buBp-on-tail unstable reduced velocity distribution has been 
c<»structed froa data obtained at the upstreao boundary of the electron 
foreshock by the GSFC electron spectrometer experiment on the lSEE-1 
satellite. This distribution is used as the initial plasma state for a 
numerical integration of the 1-D Vlasov-Maxwell system of equations. The 
integration is carried through the growth of the instability, beyond its 
saturation, and well into the stablized plasma regime. A power spectriai 
for the electric field of the stabilized plasra is computed. The spectrum 
is dominated by a narrow peak at the Bohm-Gross frequency of the unstable 
field mode but it also contains significant power at the harmonics of the 
Bohm-Gross frequency. The harmonic power is in sharp peaks which are split 
into closely spaced doublets. The fundamental peak at the Bohm-Gross 
frequency is also split, in this case into a closely spaced triplet. The 
splitting is due to slow modulations of the stabilized electric field 
oscillations which, it is thought, are caused by wave-particle trapping. 

The wave length of the m'th harmonic of the Bohm-Gross frequency is given 
by where is the wave length of the unstable mode. The mechani»n 
for excitation of the second harmonic is shown to be second order wave-wave 
coupling which takes place during that period in the evolution of the 
instability which would otherwise be called the linear-growth phase. It is 
conjectured that the higher harmonics are excited by the same mechanism. 

It is further argued that harmonic excitation at the boundary of the 
electron foreshock should be a common occurrence. 


I. Introduction 
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The theory for the generation of electrostatic plaaaa waves near the 
electron plasaa frequency upstreaa of the Earth's bow shock, as presented 
by Scarf gl.. (1971), Fredricks gjLgJL. (1971) and by Filbert and Kellogg 
(1979), has been generally accepted. This acceptance developed regardless 
of the fact that the unstable electron velocity distributions predicted by 
the theory had never been observed (Feldaun et al. . 1973) • Recently, 
observations aade on the ISEE spacecraft (Anderson gtSi.** 1981) have 
confiraed the presence of the necessary unstable distributions in the 
electron foreshock at tiaes idien significant electrostatic plasaa waves 
were being detected. A further observation of an unstable velocity 
distribution in the electron foreshock is presented here; it was obtained 
by the GSFC electron spectrometer on the ISEE-1 spacecraft (Ogilvie et al.. 
1978). Siaultaneous measurements of the electrostatic wave intensity 
obtained by the ISEE plasma wave investigation ( Gurnet t gi.. , 1978) and 

the ISEE electron density experiment (Harvey et al . . 1978) are also 
presented. These measurements indicate that the unstable distribution was 
observed just as the spacecraft was passing through the upstream boundary 
of the electron foreshock. 

The ISEE electron and plasaa wave observations both motivate and enable 
a more detailed study of the unstable plasma evolution in the electron 
foreshock. The results of such a study are presented here. 

A numerical code developed by Klimas (1982) has been used to integrate 
the Vlasov-Maxwell equations for a one-dimensional electron plasma forward 
in time from an asstaned initial plasma state. The unstable velocity 
distribution observed by the GSFC electron spectrometer was used for the 
initial velocity distribution. The initial electric field was simply 
"seeded" with very low amplitude wave modes. The results which are 
presented show the evolution of the plasma from its initial "bump-on-tail" 
unstable state through saturation of the instability and considerably 
beyond. 
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An unexpected feature of the eliMt'^ie field evolution has been 
discovered. It was expected that the saturated electric field would simply 
oscillate at the BohoMiross frequency (essent'.ally the electron plasma 
frequency) provided by the linearized electrostatic plasma dispersion 
equation (Bohm and Gross, 19^9) • It was found that not only the Bohm-Gross 
frequency but also all of its harmonics that could be included in the 
nmerical code, given its present size, were excited to significant levels 
in the saturated field. The resulting power spectria for the saturated 
electric field which is presented here shows a dMiinant peak at essentially 
the plasma frequency plus peaks near the 2fp, 3fp and 4f^ positions which 
are further split into closely spaced doublets. The splitting is due to 
low frequency modulations of the harmonic field coi^)onents presumably 
caused by wave-particle trapping effects. The peak at the plasma frequency 
is also split, in this case into three closely spaced peaks, but the 
central pesdc is so dominant that the other two are difficult to pick out. 

An explanation of the harmonic excitation is given in section IV. A 
simple verbalization of that explanation is as follows: The initial plasma 

evolution is qualitatively as predicted by linearized plasma theory (Krall 
and Trivelpiece, 1973). Field modes whose phase velocities lie on parts of 
the initial electn n velocity distribution with positive velocity and slope 
are unstable and grow exponentially with time. Other field modes Landau 
damp and therefore decay exponentially with time. All modes remain small, 
but the growing modes become very large compared to the decaying modes. At 
some point in the plasma evolution wave-particle coupling terms in the 
governing Vlasov equation which are quadratic in the growing modes become 
larger than other similar terms which are linear in the decaying modes. At 
that point the linearized description of the plasma fails and non-linear 
phenomena take over. Field modes which had earlier been decaying then are 
pumped by the dominant unstable mode which continues to grow. These piamped 
modes experience sudden shifts in their oscillation frequencies to 
harmonics of the unstable mode frequency. Interestingly, all of this 
occurs long before the saturation of the instability, during its "linear 
growth" phase when non-linear phenomena are not expected. By that time 
when the instability does saturate the pumped modes have grown considerably 
and the field is left with significant power at the harmonic frequencies. 
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The plasn aodel «ihioh ia Integrated numerioally ia a generalization of 
the uaual Vlaaov-Polaaon ayatem of equationa obtained in the eleotroatatlo 
lialt. Itezuell'a equation for the diaplacement current ia added to the 
Vlaaov (oolliaionleas Boltanann) and Polaaon equationa to provide a 
ooeplete deacrlption of the p:.aama. The generalization and the motivation 
for it are dlacuaaed in detail in aection II. Section III containa a 
diacuaaion of varioua aapecta of the moierical integration. Firat, the 
initial electron diatribution obtained from the ISEE-1 electron 
apeotroaeter experiment ia introduced. Through a compariaon of electron 
heat flux meaaurementa and aiaultaneoua plaama wave obaervatlona it la 
argued that the bxaip-on-tail unstable initial distribution was obtained 
Juat as the ISEE-1 satellite was passing through the upstream boundary of 
the electron foreshock. Next some details of the calculated electric field 
are presented and finally the electric field power spectrum discussed above 
is presented. In section IV a calculation of the second harmonic 
excitation is given which is based on a linear plasma theory with second 
order wave-wave coupling included where necessary. A comparison of the 
results of this calculation with the numerical results Indicates that the 
second harmonic excitation is due entirely to second order wave-wave 
coupling which takes place during the otherwise linear growth phase of the 
instability. It is conjectured that the higher harmonics are excited in a 
similar manner. 

II. ihe Plaaaa Model 

The results which are presented in this paper have been obtained using 
the following one-dimensional electron plasma model for the reduced 
electron distribution function, F(x,v,t), and the electric field, £(x, t) : 
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These equations are dimensionless. In the following let the starred 

quantities be the original dimensional ones. Then x s x /L where L is an 

§ 

arbitrary length scale. Solutions which are periodic in x over the 

t 

interval -L < x < L will be considered. Therefore 2L is the longest 
_ - * » 
wavelength to be considered. Time is measured by t % w^t = where 

fp is the electron plasma frequency. Then v = v /UpL and l-l /Mwen^L 

where n is the electron density averaged over the interval of periodicity; 
o 

it is a property of periodic solutions of equations 1-3 that n^ is a 
constant in t. The longest time scale of Interest is approximately 10 ms 
and the longest length scale is several tens to perhaps approximately 100 
Debye lengths. Thus, the ion density is assumed constant and equal to n^; 
the one in equation 2 stands for the dimensionless ion density. The 
magnetic fi ^Id is also assumed constant in t and only weakly dependent on 
position. The variable, x, measures position along the local magnetic 
field direction and the electric field is assumed linearly polarized in the 
magnetic field direction. The evidence that is available is consistent 
with the assumption that the electron plasma oscillations in the foreshock 
are indeed polarized along the local magnetic field (Anderson, et al . . 
1981). The weak spatial dependence of the magnetic field is ignored in 
equation 1; this is equivalent to ignoring the electron drifts. Equation 3 
is derived from Maxwell's equation, 


t+TT J + ^ X 

If V X 5 were exactly zero then U, in equation 3, would be the 

dimensionless ion velocity component along the magnetic field. Because of 

♦ ♦ 

the factor, c, in this equation, however, even small values of V x B might 
be expected to contribute. Thus, in equation 3* U is the ion velocity 
component modified by the component of c (V x B) along B. In any case U is 
treated as a constant. The net result is a model for high frequency and 
small wavelength electron plasma phenomena superimposed on a neutralizing 
ion-magnetic field background which varies only very weakly on those small 
scales. Solutions are sought subject to initial data F(x,v,0) with 
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and ^ (0) where ^ (x) is the spatial average (or zero'th Fourier aode) of 

o o 

the electric field. 

It is perhaps acre traditional when studying electron plasaa phenoaena 
to study the electrostatic lifj^t governed by the Vlasov-Poisson systea of 
equations, equations 1 and 2. But equations 1 and 2 do not deteraine 

Within the Vlasov-Poisson tradition £^(t) is normally set to zero 
and ignored. It was felt that this assioiption would be a poor one in the 
foreshock which is strongly influenced by the bow shock and hardly contains 
hoaogeneous plasaa. It would certainly be a completely untenable 
assumption when, as intended, non-periodic solutions are considered. Thus, 
the additional equation 3 has been included to ooo 4 >lete th» description of 
the plasma, including its space-averaged electric field mode. Any solution 
of equations 1-3 is also a solution of the Vlasov-Poisson system, but, with 
its space-averaged electric field governed by another of Maxwell’s 
equations. Equations 1-3 will be referred to as the Vlasov-Maxwell system 
of equa\.ions for a 1-D electron plasma. 

Klimas and Cooper (1182) have shown that any periodic solution of the 
Vlasov-Maxwell system of equations of the type discussed above can be 
transformed into a similarly periodic solution of the Vlasov-Poisson system 
with t (x) s 0 for all x, and visa-versa. The numerical integration 
results which are presented in section III are for solutions of equations 
1-3 which are appropriate, as discussed above, for the foreshock region. 

Ir section IV, however, where an analytical analysis of the numerical 
results is given, the Vlasov-Poisson system (equations 1 and 2 with £q(t) * 
0) will be used in order to take advantage of the many techniques which 
have been developed for analysis of that system. Then the transformation 
of Klimas and Cooper will be used to transform the results of that analysis 
to predictions for the behavior of the moierloal solutions. 
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111. Humerloal Reaulta 

In this section a solution of equations 1-3 will be presented using an 
initial dist ibution, derived from data obtained by the ISEE-1 Electron 
Spectrometer Experiment, which is unstable to the bump-on-tail instability. 
The numerical integration method developed by Kllmas (1982) has been used 
to obtain these results. This is a periodic solution vdiich is presented 
here as only a first step in a program of research leading to non-perlodlo 
solutions which can include the 1: fluence of the nearby bow shook in a more 
cogent manner. 

a). Electric Field 

The msnerical predictions for the electric field will be presented 
here in terns of E(t) = (L/Xp) ^(t). This scaling for the electric field 
leads to the following useful rule: When E(t) s 1, then the energy density 

in the electric field equals the total kinetic energy density in the 
Initial state of &he plasma. Since the kinetic energy of the plasma does 
not vary signifi'.antly it is true that x..en E(t) = 1 then the field and 
particle energy densities are essentially equal. In the foreshock, if this 
plasma state were reached then the electric field intensity would be 
roughly 1 V/m. 

The output of the nimierical integration scheme for the electric field 
can be writen, 

=• E-oIy) 

P 4). 

t /L-L Ay-M Cos WhITU +U'V + rttirU + UV i J 

W « I 

where 

‘ ^ + iuolol-U') Sim y] 5). 


ORIGINAL PAGE fs 
OF POOR QUALITY 


9 


in whloi) is the space average of 


uU.v') - 

The quantity, #(t), is given by, 


pOO 

\dv V \ 

— oo 



0 


The quantity, ^^(0), is part of the initial data; it oould be set to zero. 
If F(x,v,0) is such that u^(0) - U * 0, and ifpAO) is set to zero, then 
^(t) *e ^<t) a 0 and the expression for the electric field given by 
equation 4 reduces to the electrostatic limit* The solution presented here 
follows from an electron speccrometer observation for which u^(0) • U ^ 0. 
Thus, the electrostatic limit cannot be reached for this solution even if 
it is assumed that ^(0) s 0. The and are obtained directly 

from the numerical calculation. These coefficients will be presented here; 
the total field can be reconstructed using equation 4, Only ^ne field 
modes have been included in the solution which will be presented. Of 
course more modes are desirable and it is expected that more will be 
included in the future. On the other hand various solutions have been 
computed using varying numbers of modes and varying related phase 
velocities to ensure that the interpretation of the results which will be 
presented is not an artifact of the small niraber. It is i.ot expected that 
an increased number of modes will lead to any significant cLange. 


b). The Initial Data 


The dots in Figure 1 represent a reduced distribution function 
which has been constructed from GSFC electron spectrometer data taken 
aboard the ISEE>1 spacecraft on November 6, 1977 at 11:38:13-16 UT. The 
oonstructlon of the reduced distribution function requires integration over 
the velocity components perpendicular to the local magnetic field. The 
spread in the dotted curves, where it occurs, is due to the difference in 
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integration results using two methods of integration. In principle the 
integration should be carried to infinite perpendicular velocity but, of 
course, the electron spectrometer has a finite upper energy limit. The 

0 

Upper dotted ourve was obtained by extrapolating the electron spectroaeter 
data smoothly to zero for very .large velocity and the lower ourve by 
setting the electron distribution to zero for all high velocities where it 
was not measured; in each case the integration over the perpendicular 
velocity components was then carried out. Prom the results it appears that 
at least in this case, the reduced distribution is very well determined by 
the electron spectrometer data over the parallel velocity interval of 
significance to the evolution of plasma wave phenomena; i.e., the velocity 
interval containing the bump on the tail of the reduced velocity 
distribution. 

The data in Figure 1 were appF.rently collected as the spacecraft was 
passing through the boundary of the electron foreshock into interplanetary 
space. Figure 2 shows the electron heat flux determined from the electron 
spectrometer data over a one-hour interval which contains the time of the 
observation in Figure 1. The upper panel shows the heat flux magnitude, 
the middle panel the component of the heat flux vector in the direction 
measured in the ecliptic plane from the Earth-Sun direction, and the bottom 

panjl the component in the e-direction measured from the ecliptic plane. ! 

I 

The horizontal dashed lines in thr middle and bottom panels show the 

directions parallel or anti-parallel to the nominal Parker spiral magnetic 

field. The heavy horizontal line segments give one minute averaged 

magnetic field data from the ISEE data pool tapes. The heat flux vector 

can be seen to be either parallel or anti-parallel to the magnetic field f 

with abrupt transitions between. Figure 3 shows electric field data 

obtained by the Iowa Plasma Wave Experiment aboard ISEE-1 during the same { 

i 

time interval. The electron plasma frequency lies in the 31.1 kHz channel { 

during this time. Figures 2 and 3 together show the high correlation 4 

I 

between intense plasma wave noise in the vicinity of t'ie plasma frequency | 

and reversals of the heat flux vector away from its normal interplanetary 
direction alorrg the field and away from the sun. This correlation has been 
known for some time (Ogilvie gi. ftJL* i 1971; Scarf gl aJL- , 1971; Fredricks gi. 
gi. , 1971; Feldman gi. gl. , 1973) and is generally interpreted as entry into 
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the eleotrcn foreshock. The reduced distribution in Figure 1 was 

constructed from data obtain«?d during the heat flux reversal which can be 

seen in Figure 2 at apprv»xiiDatel/ 11:38. Probably because this short 

reversal occurs Just at the edge of a gap in the plasma wave date, it is 

not clearly evident in that data. Figure M contains a more detailed 

spectrogram obtained by the ISEE>1 electron density experiment over a 

fifteen minute period starting at 11;30 nnd containing the time of the 

electron observation. Tne frequency scale of the spectrogram ranges from 

essentially 0 to ^0 UHz and the grey scale indicated on the ici‘t of the 

-1 - 1/2 

figure is adjusted to cover from 16 to 72 dB above the 1 nV m Hz 
level. The short black bar on the top of the figure indicates a ti«.«. 
period during which the electron density sounder transmitter was on; this 
xs the cause of the coincident data gap in Figure 3 in the plasma wave 
data. A scan of the heat flux data and the plasma wave data in Figures 2 
and 3 shows that the satellite was in the electron foreshock at the start 
of the spectrogram. Tho general increase in the level of noise in the 
plama wave experiment 31.1 kHz channel at about 11:3^ Is clearly evident 
in the spectrogram, and in all of Figures 2, 3 and 4 the exit from the 
foreshock between 11:42 and 11:43 is very clear. The short heat flux 
reversal which starts, according to the electron spectrometer data, during 
11:38:13-16 is also clearly evident in the spectrogi'oio. The spectrogram 
contains 128 frequency steps, each of 400 Hz bandwidth, which are swept in 
a time period of 16 seconds starting at the lowest frequency and ending at 
the highest. The arrows labeled start and stop indicate the positions in 
two consecutive frequency sweeps at which the electron spectrometer 
observation of the unstable velocity distributicn began and ended. About 
six seconds before the electron spectrometer observation the satellite was 
clearly imbedded in the foreshock. About six seconds after the electron 
spectrometer observation the spectrogram indicates a quist field exct.pt for 
an intense peak at the plasma frequency. Presiasably, due to velocity 
dispersion in the foreshock boundary, by the time at which this intense 
peak is observed the satellite has passed out of the foreshock as defined 
by the low energy electrons detectable by the electron spectrometer but it 
remains in the foreshock boundary defined by the beam of higher energy 
electrons arriving from the bow shock. This beam of higher energy 
particles cont.. aes to excite plasma noise at the plasma frequency. Thus, 
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it is concluded that the lSEE-1 satellite passed out of the foreshook 
shortly after 11:38| as the data for Figure 1 were being collected, and 
soon after passed back into th^ foreshook. It is assumed that this passage 
is not clear in the plasma wave experiment data because of the da* ^ap. 

The unstable buop-on-tiil reduced distribution shown in Figure 1 is an 
example of a phenomenon which has been presumed to exist at the boundtry of 
the foreshock for quite some time (filbert and Kellogg, 1979). The 
electron spectrometer onboard the ISEE-1 spacecraft is the first experiment 
with high enough time and velocity resolution to yield that result. A 
search through the electron upectrcmeter da;a for other such exanqples is in 
progress; detailed results will be presented later. Time alising during 
tne electron spectrometer three second data collection Interval is still a 
possibility. However, at this point it appears that many of the foreshook 
boundary crossings have been resolved by the electron spectrometer and have 
been found to coincide with bump-on-tail unstable reduced distributions. 

The solid curvs in Figure 1 is a three Gaussian fit to the electron 
spectrometer data which has been used as the initial velocity distribution 
for the manerical solution to i*ollofc. This initxal data contains a cold, 
dense core, a high temperature halo, and a beam of approximately 150 eV 
electrons moving away from the bow shock along the local magnetic field. 

The electric field was assumed to initially con'^ain very low amplitude 
waves with phase velocities at the positions of the thin vertical lines in 
Figure 1. The wavelengths of the modes ranges from approximately 10 down 
to 2,5 Debye lengths (apprcximately 70 down to 18 meters). The modes with 
the largest phase velocity magnitudes are those with the longest 
wavelengths and the modes with phase velocities nearest the origin in 
velocity ir. Figure 1 are the ones with the shortest wavelength. Other 
solutions have been computed containing longer wavelength modes (Klimas, 
1982), whose phase velocities were outside the position of the beam 
velocity. Those modes were found to simply damp away with increasing time 
and to play no role in the evolution of the instability. From linear 
plasma theory, field modes with phase velocity at the position of the far 
right vertical line, on the rising portion of the bump, are expected to be 
unstable and to, therefore, grow exponentially with time until the 
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instability saturates due to non-linear phen<»ena. All other modes are 
expected to decay exponentially with time (to Landau damp). The aotual 
evolution of the field modes turned out to be considerably more complex. 

c). The Solution 

Figures 5a-d contain plots of the absolute values of 6^(t) through 
B^(t) over more than one hundred plasma periods (approximately 3 ns. total 
elapsed time) starting at the initial plasma state. The A-coeffloients 
differ from the B-ooefficlents only in details and therefore space will not 
be taken to present plots of them. A close examination of these figures 
will reveal that the curves are made up of a sequence of peaks of slowly 
varying amplitude. Each pair of peaks represents a single period of 
harmonic oscillation on this semi -log plot of the absolute values. Figure 
5a, for B^(t), shows a classic example of instability growth and 
saturation. Following an initial transient the anq?lltude of B^(t) grows 
exponentially at a rate which can be predicted very accurately using a 
dispersion relation fr<XB linear plasma theory. The oscillation frequency 
of this mode is also very accurately predicted by the Bohm-Gross relation 
from linear plasma theory. (These issues will be discussed in detail in 
the next section.) The instability saturates at t > 200 leaving this mode 
oscillating at the Bohm-Gross frequency, very close to the plasma 
frequency, with slow modulations of the amplitude which are due presumably 
to particle trapping effects. 

A detailed description of the computed electric field will be given in 
the following subsection. At this point, however, a brief argument can be 
given which shows that the dominant feature of the electric field is in 
good agreement with the field detected by the plasma wave experiment: The 

total field, following saturation, is dominated in the vicinity of the 
plasma frequency by the contributions of B^(t) and A^(t). From this 
mmerical solution it would be expected that plasma wave noise near the 
plasma frequency witn anvlltude “ 100 mV/m and wavelength close to 70 
meters should be detected by the plasma wave experiment. The calibration 
of the :las»a wave experiment is such that the pl< t of the 31.1 kHz channel 
in Figure 3, which contains the plasma frequency, should reach full scale 
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•t several tens of mlllivolts/Beter (private oonnunioation, D. k. Qumett, 
P.I., Plasma Have Experiment)* An examination of Figure 3 reveals that the 
period Just preceding and following the time of the electron spectrometer 
observation is oh*>**aoterized by unusually strong plasma wave activity near 
the plasma frequency. The 31*1 kHz channel appears full scale or near full 
scale almost continuously. The 70 meter wavelength predicted by the 
numerical solution is smaller by about a factor of three than the 215 meter 
antenna used to obtain the data in Figure 3. Thus, the plasma wave 
experiment would be expected to detect an electric field oscillating near 
the plasma frequency with amplitude roughly 30 mV/m, Just about its full 
scale field level snd in agreement with the data presented in Figure 3. 


The plot of BgCt) in Figure 5b provides the first evidence of quite 
unexpected evolution for this bump-on-tail instability. Initially BjCt) 
decays exponentially in time as expected with rate close to the Landau 
damping rate predicted by linear plasma theory. The reversal of this decay 
at T > 70, the subsequent exponential growth, and the simultaneous shift to 
a higher oscillation frequency were not expected. After all, none of the 
usual assvunptions which one maxes to Justify the accuracy of the linear 
plasma theory appear violated during this reversal but it is clear that the 
linear plasma theory (cr even quasi-linear theory (Davidson, 1972)) would 
not allow such behavior. This situation is made even more difficult by the 
plots of B^Ct) and B^(t) in Figures 5c and 5d which show more dramatic 
examples of the same kind of behavior. It would appear that the number of 
modes involved In this behavior is limited here by the number of modes 
included in the calculation and not by any feature of the physical 
phenomenon involved. On the other hand it is important to realize that not 
all possible field modes are involved. Other numerical solutions have been 
computed containing electric field modes with wavelengths longer than that 
of the unstable mode in this run; those field modes did not experience a 
reversal of their Landau damping. In addition, it has been found that 
field modes with wavelengths shorter than the unstable mode which, however, 
do not satisfy = x /m, where X is the unstable mode wavelength and m is 
any integer greater than one, also do not take part in the decay reversal 
and frequency shift. The solution presented here has been arranged so that 
t.ie unstable mode .’s at the longest wavelength, all other modes are at 
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shorter wavelengths, and they all satisfy s x„/b. Thus, in this 

B U 

solution, only those modes that can reverse have been Included because. In 
addition to the unstable mode, those are the only modes which rise to high 
enough amplitudes to make measurable contributions to the electric field. 

d). The Electric Field Power Spectrum 

The spectral distribution of the computed electric field is 
considerably richer than was expected. The electric field is not limited 
to plasma oscillations at (or very near) the plasma frequency. The decay 
reversals and simultaneous shifts to higher frequencies discussed above 
lead to a power spectrum with peaks at all the hamtonics (limited only by 
the numer of modes in the numerical calculation) of the fundamental 
Bohm-Gross frequency. In addition these peaks are split into doublets or 
triplets due to the aBg>litude modulations evident in Figures 5a-d following 
saturation of the plasma instability at t » 200. 

Figure 6 contains the electric field power spectrum computed from the 
output of the numerical calculation and, also, a power spectrum of a model 
electric field which is used to help interpret the numerical results. The 
electric field was computed using equation 4 with U s 0; this is the field 
in the reference frame moving with the solar wind. It will be shown bel<w 
that the wavelengths in the electric field are quite small. As a result, 
the Doppler shifts associated with the transformation from solar wind to 
satellite frame of reference turn out to be very small corrections to the 
frequencies of the important field modes and, in general, can be ignored. 
The constant, £^(0), nas been set to zero in order to present the simplest 
case. The vertical tick marks at the top of the figure give the positions 
of the plasma frequency and two, three and four times the plasma frequency. 
The spectrum was computed from that portion of the solution following the 
saturation of the instability, from t = 200 to t “ 630. The field prior to 
t “ 200 is relatively much weaker and after 200 shows no evidence for 
further qualitative evolution. Thus, it is felt that this portion of the 
solution yields the electric field that the plasma wave experiment should 
detect when it detects peaks in the electric field at the foreshock 
bounvdary. The heavy curve gives the power spectrum computed from the 
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maerioal solution. Bsoause the tiae interval over which the eleotrlo 
field was ooaputed is short there Is oonsiderable leakage in this speotrua* 
To evaluate this effect a power speotrue for a aodel of the nue>rioal field 
was also computed, k model eleotrlo field with pure harmonics at the 
Boha-Gross frequency of the vmstable mode and at two, three and four times 
that frequency was constructed. The amplitudes of each of the haraonios 
were chosen using the following reasoning: It was found that each of the 

well separated groups of peaks in the p<wer spectrin came from distinct 
wavelengths i.~i the field. The peaks in the vicinity of the Bohmmross 
frequency corresponded to the longest wavelength, or the m s 1 terms in 
equation 4. The peaks in the vicinity of twice the Bohm-Qross frequency 
come from the next shorter wavelength, half the longest, or the m « 2 teras 
in equation 4. And so on. The following general rule seems to apply. 

After the decay reversal and frequency shift at x « 70, each cf the field 
modes that satisfy s x^/m experiences a shift in oscillation frequency 
such that its phase velocity, moves up to that of the unstable mode; 

i.e., up to the rising portion of the biun? on the reduced velocity 
distribution. All of these modes then have the same phase velocity and 
their oscillation frequencies are given by = m where is the 

Bohm-Gross frequency of the unstable mode. None of these modes with m > 1 
satisfy a dispersion relation; they are not normal modes of the plamsa. 
Because each group of peaks in the power spectruD had its source in a 
single m-value in the sum of equation 4 it became a relatively simple 
matter to choose the amplitudes of the harmonics in the model electric 
field. The electric field contributions of each of the m-values in the sva 
were scanned for a maximum value. These maxima were then used for the 
amplitudes of their respective (same m>value) pure harmonics. The power 
spectrvmi of this model field made up of pure harmonic oscillations was then 
computed; the thin curve in Figure 6 is the result. The overall shift 
upward from the heavy to the thin curve is due to the choice of the maxima 
for the amplitudes of the pure harmonics. It is clear from the OKidel field 
spectrum that there is indeed considerable leakage. The true spectrum for 
the model field, if ’ , were computed over an infinite time interval, would 
consist of spikes at the Bohm-Gross frequency and its harmonics with no 
power in between. It is also clear from Figure 6 that the computed field 
power spectrum is consistent with power in isolated peaks each of which are 
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split into several closely spaced co^>onents. The second, third and fourth 
haraonics are clearly split into doublets and the first haraonic appears to 
be split into a triplet. This splitting of the peaks in the power speotrua 
is due to the slow aaplitude aodulations evident in Figures 4a-d after t • 
200 which are due probablj to particle trapping in the electric field. 

IV. A Mechanlaa for Higher Haragnic ExaiiiaUon 

It is thought that the decay reversals and frequency shifts discussed 
in the previous section can be explained in teras of a second order 
wave-wave coupling aechanisa. This mechanism allows the unstable plana 
aode to pump the stable ones to significant amplitudes with oscillation 
frequencies that are narmonics of the unstable mode frequency. A specific 
calculation of the amplitude, growth rate, and oscillation frequency of the 
second harmonic is given in this section. The results of this calculation 
are in excellent agreement with the numerical computation. It is 
conjectured that the higher harmonics are due to the same mechanism 
involving many more wave-wave couplings but a detailed calculation of the 
higher hamionics has not been done. 

The calculation of the second harmonio excitation is essentially a 
linear calculation with second order coupling added where appropriate. As 
discussed above, the calculation will be done in the electrostatic limit 
(i.e. , solutions of equations 1 and 2 with the constraint Sq(t) = 0 will be 
studied) and the results will then be transformed to allow a comparison 
with the numerical computation. Therefore, consider solutions of equations 
1 and 2 with£^(x) = 0 which are periodic in x over the interval - 1 < * ^ 

1 and are subject to Initial data F(x,v,0). As in the numerical 
computation, let F(x,v,t) and £ (x,x) be represented by truncated Fourier 
series expansions in x with coefficients given by 

A I C “ Civi Tl X _ 

= 5 y A e 1- U,v^ V) 

- I 

and 

I r 

= A \(lx £ 




5 (x , v ) 
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Then, these coefficients satisfy 

^ ■ f r Ua 

ir + LmT\vV„ * /_ JV 7). 

• n=” 

^vn ~ wnr «)• 

.. oo 

and 6^ s 0. Due to the truncation equatio.is 7 and 8 govern those 
coefficients which satisfy jsj < M and in the siounatlon in equation 7 |m-n| 
< M and |n| < M where M is some integer greater than zero. In the 
numerical results presented in the previous section and for the rest of 
this calculation M s 4. All other coefficients are assumed zero. 

To proceed assume that the plasma contains only very weak plasma waves, 
i.e., assume that f << f for m ^ 0 and assume further that the f are 
well behaved as Iv) ■*■<*> such that thet calculated from equation 8 are 
similarly small. These are standard assumptions which one would make in 
order to proceed with a linear theoretical treatment of equations 7 and 8. 
This will be the approach taken here until an obvious contradiction to the 
linear theory is found. Since = fcju only m > 0 will be considered. 
Thus, consider the linear approximation to equations 7 and 8 starting with 
m = 0. 

a) . m s 0 

r = 0 ’>• 

There are no linear terms in the summation of equation 7 because 6^ - 

0. Thus, f (v,t) = f (v,0) = g(v). 
o o 

b) . m = 1 

This is one of the unstable modes. The linear approximation to 
equation 7 is 
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Of 'lOursA the problM of finding a solution to equations 8 and 10 for f^ 
and€^ is e faailiar one; the approach taken here is due originally to 
Landau (19*l6). Let f^(v,p) and l^^(p) be the Laplace transforms of f^ and 
6^ with the symbol, p, for the Laplace variable. Then, 
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With f (v,0) 

and f^(v,0) given by the 

sum of three gaussians 

in V, as 

in 

the case here, the method of Jackson 
D.(p) which leads to the simple pole 

(I960) can be used to find the zero in 

/V 

in G^(p) that domin*ites €^(t) as t 

grows large. 

With, 


a 





— 

p-po 


14). 


where p^ s o 4- l u, the result of this calculation is o s .0481 , a positive 
growth rate and therefore an unstable mode, and u s >1.07 (remember that in 
these dimensionless units |o>| & i corresponds to oscillations at the pla«aa 
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frequency). Generally, the method of Jackson, when It applies, gives a 
more accurate result for u than does the Bohm-Gross expression. However, 
in this particular case, with m s i, the two results are essentially 
Identical and therefore this frequency has been referred to everywhere else 
in this paper as the Bohm-Gross frequency of the unstable mode. The lowest 
frequency peak in the model field power spectrum (the thin curve) of Figure 
6 is at this calculated Bohm-Gross frequency. The excellent agreement 
between the position of this peak and that derived from the nunerlcally 
computed field (the heavy curve) indicates tUe accuracy with which the 
linear theory predicts the oscillation frequency of this unstable mode. 
Figure 7 shows the details of B^(t) during its linear growth phase. The 
slope of the solid line drawn on that figure can be used to calculate the 
actual growth rate of the instability; the result is .0530 as opposed to 
the calculated o = .0481. Thus, the growth rate of the unstable mode is 
also very well rredicted by the linear theory. Generally, the behavior of 
the unstable modes during their linear growth phase is very well predicted 
by the linear plasma theory. 

c) . m = 2 

This is the first of the modes which experience the reversal of their 
Landau damping and the subsequent growth and frequency shift which are 
" ith in violation of the linear theory. It will be shown here that the 
evolution of this mode during its decay and subsequent growth is governed 
by 


ly ■** ^ ^5). 

At t s 0, fp and f^ are equal in amplitude as areS 2 C^. All of these 
quantities are very small. Thus, the last term in this equation, which is 
quadratic in these small quantities, is initially completely negligible. 
With the neglect of this last term, ^otii and C ^ obey the same linear 
theory discussed above in the m = 1 case. A similar calculation of o and u 
in this case yields o = -.151 and w = -1.27 in excellent agreement with the 
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behavior Indicated in Figure 8 of B 2 (t} during its Initial decay. However, 
the neglect of the last term in equation 15 cannot remain valid. While the 
linear theory applies and f ^ are growing exponentially in time and 
and are decaying exponentially in time. At some point the last term 
must dominate the right-hand side of this equation. In this subsection the 
effects of the last tern, when it dominates, are calculated. It is shown 
that the decay reversal and frequency shift of B 2 (t) are thereby completely 
explained. 
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Using equations 15 and 8, it can be shown that the Laplace transform of 
(t) is given by. 




L 

ATTO;,(p^ 



r v-i I f ^ ft/i ‘ 

• — oo 

in which f* = af /dv and the path of integration, C, is from -i» to 

^ ^ * iT 

with €^(p') analytic to the right and f'-j(v, p-p') analytic to the left of 
the path on the complex p'-plane. The first term in this equation would be 
the sole result if a linear calculation were being done; it predicts the 
initial behavior of BgCf) while it is decaying. Let " ^^P^ ^^P^ 

where £(p) is the linear first term in equation 16 and Q(p) is the 
quadratic second term. With equation 14 for 


'W ^ i I 0^ \ ^ 

Q'P'' = ilTD^ip^ )c*v'(pvUlTv’l S'lVjp-p.) 

' ->oe 

If equation 11 is used for then it can be shown that ($(p} contains one 
term with a pole at p a 2p^; this tarn is given by, 
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Finally, the contribution of this pole to i* given by. 
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This expression gives a contribution to 62(7) which oscillates at exactly 
twice u, the Bohm-Gross frequency of£^, and grows exponentially with twice 
the growth rate, 0, of B^. The amplitude of this contribution to BgCx) is 
completely determined by the fit, in Figure 7, to the amplitude of6 ^(t) 
which yields |a| and a and by the calculation of u done in the previous 
subsection; there are no free parameters available. The prediction from 
equation 17 for the amplitude of BgCx) is given by the straight line in 
Figure 8. The prediction from equation 17 for the oscillation frequency, 
2w, is given by the position of the second peak in the model field power 
spectrum in Figure 6 which lies Just in the middle of the corresponding 
doublet determined from the numerical solution. It would appear that this 
contribution to^2^‘^) accounts entirely for the decay reversal and 
frequency shift to the second harmonic in B2 (t) and thus, for the second 
harmonic excitation in the electric field. The mechanism for this 
contribution which has been examined above would normally be called second 
order wave-wave coupling but this is, perhaps, somewhat confusing 
terminology si<~ce it is actually the m s 1 mode interacting with itself 
which produces this result; it is hard to imagine a wave interacting with 
itself. Actually, the interaction is between the m s 1 wave in the 
electric field and the m = 1 density fluctuations in the particles; the 
terminology is simply inappropriate. 


It is conjectured that the higher harmonics in the electric field are 
produced by the same mechanism as that which produces the second harmonic 
but, with more than Just one mode doing the pumping. For example, after 
its decay reversal the third mode shifts to the third harmonic of the 
fundamental Bohm-Gross frequency and grows exponentially with a growth rate 
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which is three times that of the first mode. It is conjectured that this 
mode is governed by equation 8 and 




v.3ffv5j 




The first term on the right-hand side of this equation would provide the 
initial linear decay of £^(t) and the last two would do the pumping, k, 
the index, m, increases the number of possible mode-mode couplings 
responsible f-.r the pumping goes up rapidly. 


V. Conclusion 


The essratial element which allowed the excitation of the harmonics in 
the example discussed in this paper was the large ratio of free energy in 
the electron beam to the field energy prior to the growth of the 
bump-on- tail instability. Thus, the unstable electric field modes grew by 
several orders of magnitude before the free energy was exhausted thereby 
leading to the important consequence that the linear growth phase of the 
instability lasted for a long time. During this long linear growth phase 
other shorter wavelength field modes which were Landau damping were able to 
decrease their strengths by several orders of magnitude. The result was a 
large imbalance between the strengths of the stable and unstable modes in 
the plasma; so large, in fact, that the linear Landau damping of the stable 
modes was eventually dominated by wave-wave coupling quadratic in the 
unstable modes. Although a quantitative statement cannot be m> de, it does 
seem safe to say that the thermal electric field noise level in the 
undisturbed interplanetary plasma at the frequencies and wavelengths of 

-5 

Interest here is very small, probably smaller than the approximately 10 
V/m which was used to seed the field modes in the numerical computation 
(Meyer-Vernet , 1979; Hoang £La 1., 1980; Couturier fii. ai.. , 1981). Thus, 
the free energy of an electron beam at the foreshock boundary should easily 
dcmiinate the undisturbed field energy. Certainly the plasma wave 
observations at the boundary of the foreshock are often consistent with a 
growth of several orders of magnitude in the field strength during the 
evolution of the instablity. It appears then that the excitation of 
harmonics at the boundary of the electron foreshock should be a common 


occurrence. 
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The nuaerloal results which haws been presented ha\e been ooiq>uted 

using only a raall number of discrete field modes. However, various 

unreported nunerioal solutions tiave been ooq>uted to test the conclusions 

which have been given. It has been found that field modes with wavelengths 

longer than that of the unstable mode do not play any detectable role in 

the evolution of the unstable mode or the pi^ed nodes. It has also been 

found that the sane can be said of field modes with wavelengths shorter 

than that of the unstable mode whose wavelengths, however, do not satisfy 

the equation X s X /a. The solution which has been presented contains 
n u 

only field modes which do satisfy this equation with 1 < a < ^. It is 
thought that only modes which satisfy this equation can possibly have a 
significant part in the evolution of the Instability. Thus, it is thought 
that the solution which has been presented is essentially ooBq>lete except 
for the modes which satisfy x. » X /m with m > 4. But the wavelengths of 
those modes, according to this equation, would be even shorter than the 
minimum 2.5 X^ wavelength which is Included here. It is not considered 
likely that field modes with wavelengths equal to Xp or shorter would play 
a significant role if they were included but this point will be checked in 
the future when solutions with more field modes will be couq^uted. 

It is further not considered likely that the transition fr«B the 
discrete modes appropriate for the nuserioal oooq>utatlon to a continuum of 
modes more appropriate for an Interpretation of the plasma wave data would 
have any significant impact on the harmonic excitation. Following the long 
linear growth phase of the Instability the contribution to the electric 
field from the continuum of unstable wavelengths should be dominated by 
those wavelengths, which ere very close to the most unstable wavelength. 
Because the field would therefore be dominated by a narrow range of 
wavelengths, and because the wavelength dependent Bohm-Gross correction to 
the plasma frequency is small, all of the important unstable field modes 
should oscillate at very nearly the same frequency, near the plasma 
frequency. With the resulting sharply peaked wavelength and frequency 
spectra for the unstable field modes it would be expected that the spectra 
for the pumped harmonic modes would be similarly peaked. Thus, it is 
expected that in the continuum, during the linear growth phase of the 
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Instability and nhilt tha hanaonio puaping is ooourring, tha fiald would 
avolva in a aannar quits siailar to that of tha diaorata nods oaloulation. 

Than tdiy hava tha haraonios of tha plana fraquanoy not baan obaarvad 
at tha alaotron foraahook boundary? Tha anawar ny ba that thay hava baan 
obaarvad but dua to tha natura of tha plana wava dataotora it is difficult 
to dataraina that this is tha oaaa. An unfortunata aids affaot of tha high 
sanaitivity of thasa dataotora is thair tandanoy to spill strong signals in 
any ona fraquanoy ohannal ovar into othar naarby ohannals. A good axnpla 
of this phanoaianon oan ba saan in Figura 3 in which a very strong signal at 
tha plaan fraquanoy in tha 31*1 kHz ohannal ny ba axoiting uny of tha 
othar fraquanoy ohannals in tha dataotor. Thara ny also ba harnonios of 
tha plana fraquanoy in that data but that would ba vary difficult to 
prova. This is an axtram axnpla; thara ny ba othar tinea tdian harnonios 
in tha pljaaa waves oould ba separated fron spillover in the detector. On 
the basis of the results of this paper it would appear that thay are thara 
to be found. 
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Figure 1 The dots indicate an electron reduced velocity distribution 
constructed froa data obtained by the electron speotrooeter 
experiment on ISEE-1 as it passed through the boundary of the 
electron foreshock. The solid curve is a three gaussian fit 
used for initial data for the numerical integration. The 
thin vertical lines indicate the initial phase velocities of 
the field modes included in the numerical integration. 


Figure 2 Electron heat flux determined through a moment fit to data 
from the electron spectrometer experiment on ISEE-1 over a 
one hour ..nterval which contains the time of the observation 
in Figure 1. The sudden shifts by approximately 180* in the 
angle ^ indicate passage through the boundary of the electron 
foreshock. 



Figure 3 Electric field measurements made by the plasma wave 

experiment on ISEE-1 using the 215 m. antenna over the same 
time interval contained in Figure 2. The plasma frequency 
during this time interval is in the 31.1 kHz channel. Sudden 
increases in the signal in this channel can be seen to be 
very well cor.'elat d with the heat flux reversals in Figure 
2. Tlid periodic data gaps occur when the electron density 
sounder experiment is turned on. 

Figure U Electric field spectrogram made from data obtained by the 

electron density experiment on ISEE-1 over a fifteen minute 
time interval containing the time of the electron 
spectrometer observation in Figure 1. The spectrogram is 
constructed frcxn consecutive 16 second frequency sweeps 
starting at the lowest frequency and ending at the highest. 
The beginning and end of the electron spectrometer 
observation are at the points indicated by the arrows labeled 
start and stop. 
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Figure 5 Four of the electric field modes obtained from the numerical 
integration of the 1-D Vlasov-Naxwell equations using the 
bump-on-tail unstable reduced velocity distribution given in 
Figure 1 for the initial plasma state. The mode labeled 
contains the unstable fundamental. The other three modes are 
initially stable but are later pimped by the uni;table 
fundamental. The pumped modes oscillate at frequencies which 
are very close to the harmonics of the plasma frequency. 

Figure 6 The heavy curve is a power spectrum constructed from the 
computed electric field following saturation of the 
bump-on-cail instability. The light curve is a power 
spectrum constructed in a similar manner from a model 
electric field containing power only at the Bohm-Gross 
frequency of the fundamental unstable mode and its harmonics; 
the power apparent in this curve between the peaks is due 
solely to leakage. A comparison Indicates that the computed 
electric field contains power in isolated peaks, at the 
harmonics of the Bohm-Gross frequency, which are further 
split into closely spaced components. 

Figure 7 The linear growth phase of the mode plotted in Figure 5a. 

The straight line is a fit which has been made to determine 
the magnitude of the oscillations in this mode and their 
growth rate. The growth rate determined in this manner is in 
very close agreement to the growth rate predicted by linear 
plasma theory. 


Figure 8 The B^ mode from Figure 5b during the linear growth phase of 
the burop-on-tall instability. The initial decay is as 
expected from linear plasma theory. The reversal of that 
decay and simultaneous frequency shift are due to pimping by 
the fundamental unstable mode through quadratic wave-wave 
coupling. The straight line is a no-free-parameter fit to 
the amplitude of this mode, during the pumping phase, 
calculated using the wave-wave coupling mechanism. 
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